Reconstruction of the electron-phonon interaction function in Ta, 2H-NbSe 2 and 
MgB 2 from spectra of S — c — N point-contacts 



N. L. Bobrov *, V. V. Fisun, 0. E. Kvitnitskaya, V. N. Chernobay, and I. K. Yanson 

B.I. Verkin Institute for Low Temperature Physics and Engineering, 
of the National Academy of Sciences of Ukraine, prospekt Lenina, ^7, Kharkov 61103, Ukraine 
(Dated: July 30, 2012; Published Fiz. Nizk. Femp. 2012, 38, p.480(Low Femp. Phys. 2012, 38, p.373)) 

Using Fa, 2H-NbSe2 and MgP>2 as an example it is shown that it is possible to reconstruct qualita- 
tively a function of the electron-phonon interaction from point-contact spectra in a superconducting 
state. The limits and the restrictions of this method are also shown. The results obtained are 
compared with available literature data. 
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I. INTRODUCTION 

A traditional point-contact (PC) spectroscopy of the 
electron-phonon interaction (EPI) 1 assumes that both 
electrodes forming a N — c — iV-contact are in a normal 
state at low temperatures. However, it is quite difficult 
to transform many actual superconductors to a normal 
state. Spectra of S — c—N and S— c— S'-point-contacts in 
some cases contain an additional structure in the region 
of phonon energies, and a form of such spectra is quite 
complicated and on occasion is far from traditional. A 
transformation into a normal state of electrodes form- 
ing a point-contact often leads to disappearance of such 
a structure. It can be explained by the fact that the 
spectrum in a normal state and the additional nonlin- 
earity arising at a transition to a superconducting state 
are determined by scattering processes in geometrically 
different regions of a point-contact. 

In a normal state the contribution to a spectrum is 
backs cattering processes at which an electron returns to 
the same half space from which it flew oulji. A main con- 
tribution to a spectrum is formed by processes of scatter- 
ing of electrons in a volume with a diameter of the order 
of a point-contact. Indeed, for a contact in the model 
of an aperture with a diameter d a probability of the re- 
turn in the case of the scattering is maximal in a center 
of the contact and is 0.5. Since in an isotropic case a 
direction of the electron flight after scattering is equally 
possible, on a large distance R ^> d from the contact 
plane the probability of the return does not exceed a ra- 
tio of the contact area to the surface area of a sphere with 
a center at the scattering point and the radius R, and is 
^(d/4R) 2 . Thus, on the distance R ~ 2c? this probabil- 
ity is 32 times lower than in the contact center, and 
for R ~ Ad it is ~ 128 times lower. A maximal con- 
centration of nonequilibrium phonons is achieved in the 
region of maximal density of the current and decreases 
fast with its spreading. So in the first approximation at 
large distances from a narrowing the point-contact can 
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be considered as a point source of phonons, and their 
density decreases ~ 1/R 2 . However, the largest part of 
phonons is still generated in edges in the amount of the 
order of the energy relaxation length of nonequilibrium 
electrons, which depends on their energy. The inelas- 
tic mean free path of electrons can be estimated by the 
formula l £ = vft £ , 

here g{u>) is the Eliashberg thermodynamic function. 
Note that if a metal has a few groups of electrons, then 
in the estimation it is necessary to use the Fermi velocity 
for the group which determines its transport character- 
istics. By the way, this note is also valid for the expres- 
sions given in Sec. II for an estimation of the potential 
barrier appearing at an interface of two metals due to 
a mismatch of their Fermi parameters. Obviously, that 
the scattering of nonequilibrium electrons by phonons in 
edges docs not affect the PC spectrum in a normal state. 
An inelastic superconducting addition in the spectrum is 
determined by processes of scattering of electrons, under- 
going the Andreev reflection, by nonequilibrium phonons, 
and any scattering process is effective. Since the proba- 
bility of reabsorption of nonequilibrium phonons by elec- 
trons experiencing the Andreev reflection depends on the 
concentration of phonons, from the above estimation it is 
clear that the contribution to the superconducting inelas- 
tic spectrum decays more slowly with increasing the dis- 
tance from the aperture than the contribution to the nor- 
mal spectrum. Besides, not taking into account phonons 
created in edges would be too rough approximation, since 
the Andreev reflection takes place in a volume of the or- 
der of the coherence length, which is, in the cases under 
discussion, considerably larger than the diameter of the 
point-contact, and is comparable to the inelastic relax- 
ation length at the Debye energies. Estimations for the 
concentration of nonequilibrium phonons depending on 
an applied voltage are reported in Ref|2j for a Nb point- 
contact. It seems that a relative integral contribution of 
the peripheral regions of a contact in the inelastic super- 
conducting addition in the spectrum for superconductors 
with a sufficiently strong EPI can be comparable with 
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the contribution of the central regions. Therefore, if the 
contribution to the spectrum of the point-contact central 
regions is "turned off' somehow, the remained nonlinear- 
ity will almost entirely be determined by the supercon- 
ducting addition to the spectrum of peripheral regions. 
Taking into account that the impurities and defects dis- 
turbing the crystal structure of the material under inves- 
tigation are usually concentrated near the surface, and 
considering that in the process of forming point-contacts 
additional distortions of the lattice are introduced, sig- 
nificantly lower quality of the material near the contact 
center can be expected. This leads to a blurring of the 
normal spectrum and a decrease in its intensity, and to 
the complete suppression of it for a very strong distortion 
of the lattice. At the same time edges of the contact can 
remain quite clean, and the superconducting addition to 
the spectrum will reflect properties of the bulk material. 
We note in passing that the observed high level of the 
background for some PC spectra in a normal state, which 
cannot be described within the framework of the exist- 
ing theory of reabsorption of nonequilibrium phonons 3 , 
possibly is also related to the presence of a strongly dis- 
torted layer with a thickness smaller than the diameter 
of the contact. The aim of this work is reconstruction 
of the PC function of the EPI from the superconducting 
addition to the spectrum using clean point-contacts of 
Ta as an example. Then, using algorithms obtained for 
Ta, the functions of the EPI for 2H-NbSe 2 and MgB 2 are 
restored, spectra of which in a normal state are absent 
or blurred very much. 



result coincides with that obtained in RefH 
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Taking into account an effect of the EPI on a form of the 
IVC at high eV gives nonlinear in V contributions to the 
current, 

I(V) = I + 6I&(V) + I° exc + 5I* h (V) 

in this case, 8I^ h (V) coincides with a negative addition 
in the current obtained for a normal contact 1 , 



8I£ h (V) = -4eN(0)n / duL(u,eV,T)G"(u) (5) 



Here fi ~ d 3 is the volume of generation of phonons, 



L(oj,eV,T) = -(eV-oj) 



cth 



eV - oj 
2T 



cth 



CO 

"2T 



(6) 



The superconducting component of the inelastic current, 

8I$ h (y) = -4eQ7V(0) [f^°dojL 1 (oj,eV,T,A)G N (oj) + 
+ / °° dojL 2 (oj, eV, T, A)G S M] 

(7) 

For eV - oj^> max(A, T) 



Li(w) 



8A 



L 2 (oj) 



2A 

IT 



(8) 



II. MAIN THEORETICAL CONCEPTIONS 

The authors of Ref|4| have extended a theory of the PC 
spectroscopy to S — c — S contacts. In this case it has 
been assumed that S — c—S contacts satisfy the following 
conditions: 
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These conditions mean that characteristic dimensions 
of the contact are smaller than the momentum mean free 
path li and the coherence length £(0), the energy mean 
free path l £ 3> d up to energies Tiuid, (l E ~ vf/^d for 
e ~ fiuJn)- The I — V characteristics (IVC) in the zero 
approximation in <f/£, d/l e is given by the expression, 
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where R^ 1 = e 2 p 2 F S /it 2 is the resistance in the N state 
for a quadratic dispersion law; S is the aperture area. In 
the case of a real value of A for an excess current the 



The function G N {u) in Eqs.© and (J7J is a PC function 
of the EPI which is similar to that in point-contacts in a 
normal state, G s (ui) is a superconducting PC function of 
the EPI different from G N (w) by a form- factor. In this 
case, in contrast to a normal form-factor causing the con- 
tribution to the current due to electron-phonon collisions 
with changing the z-component of electrons velocity, in 
the case of the superconducting form-factor entering into 
G s (oj) a contribution to the current is made by electron- 
phonon collisions related to processes of the type of the 
Andreev reflection in the contact region, i.e. the transfor- 
mation of quasielectron excitations into quasihole ones. 
A relative value of the phonon contribution to the ex- 
cess current at eV ~ ujb is of the order of d ■ ujf>/vf, i.e. 
it is small if the condition <j2j) is fulfilled. For the first 
derivative of the excess current it is obtained the relation 
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Since as shown in Ref. at eV -C ojd, G n (oj) ~ 
(uj/uj£>) 4 , G s ~ oj, then in the region of small displace- 
ments on the contact a main contribution is made by 
processes without changing the z-component of the ve- 
locity, of the type of the Andreev reflection. In Ref. 6 
S — c — N point-contacts have been considered. For the 
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second derivative of the IVC it was obtained the following 
expression 
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where G a {ui) is the function of the EPI for normal and 
superconducting metals forming a heterocontact, S(x) is 
the blurring factor, 



S(x) = 9(x - 1) 
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Thus, at T — > a resolution is determined by A. For 
the first derivative of the excess current, the expression 
similar to the S — c — S contact was obtained, 
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As it was mentioned above the expression for the excess 
current of the pure S — c — S contact, similar to (jlj, was 
obtained in Refi. In that work as well as in RefQ, a 
form of the IVC for S — c — S and S — c — N contacts 
in pure and dirty limits was found. In a dirty limit the 
excess current is lower and at arbitrary temperature for 
the S — c — S contact is described by the expression^ 
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where U is the momentum mean free path, d is the con- 
tact diameter. In the intermediate case, for an arbitrary 
relation between Zj and d an evaluated relation, correct to 
first order in d/U, was obtained in RefSfor the S~c-S 
contact, 
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Here R^ 1 = R^ 1 (1 — d/li) is the contact resistance 
in a normal state with taking into account the impurity 
scattering, Rq 1 = e 2 SSF (a) (27r?i) -3 - is the resistance 
of an extremely pure contact. An absolute value of the 
excess current for S — c — N contacts appears to be twice 
smaller. In Refs. 0-H upon calculating the excess cur- 
rent, it has been assumed that the Fermi velocities and 
momenta coincide at both edges of the point-contact. An 
effect of a barrier within the plane of narrowing, appear- 
ing for different vf and pf, has been taken into account 
in RefH. For the excess current of a pure S — c — N het- 
erocontact at arbitrary temperature it was obtained the 
expression, 



j(0) 



Rn 



Jth(V/2T) 



(15) 



J = 



2(aD) 



D 2 



i - 



D 2 arcthVl 
71' (HTR) 



(16) 



Here R — 1 — D is the reflection coefficient, D = D{a) 
is the diffusion coefficient (calculation of D(a) and other 
coefficients for heterocontacts is considered in detail in 
Ref. H). For fl->U = 4/3. With decreasing the diffu- 
sion coefficient the form of the IVC and the dependence 
of the differential conductivity S — dl/dV on V approach 
dependences typical for tunneling contacts. The reflec- 
tion from an interface between metals leads to that S(V) 
at low temperatures becomes to be non-monotonic. At 
T = 0, 



(aD(a)) = a N /a(A) 
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Thus the transparency can be estimated from experimen- 
tal data. For the S — c— S contact for an arbitrary trans- 
parency of the interface in the region of high voltages 
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where J is determined by Eq. (jl6|) . As it was already 
mentioned above, the second derivative of the IVC for 
the S — c — N point-contact is given by the convolution 
Eq. (jTU]) . Taking into account the relation between the 
derivative of the IVC and the PC function of the EPI one 
can write down 
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As a model gp C (to) we will use a function of the EPI 
for tantalum, restored from a spectrum of a heterocontact 
in a normal state. Note that on writing the formula (|T9]) 
we neglected the fact that G a (to), corresponding to the 
superconducting edge, is a superposition of g N (uj) and 
G s (ui) (see Eq.©). As it was mentioned above, G s (uj) 
is different from g N (ui) by a form-factor, and most pro- 
nounced differences are expected in the region of small 
displacements. In what follows, since we obtained the 
spectrum of the EPI in a superconducting state for the 
same contact, a comparison of theory and experiment 
will be performed, and it will be considered how much 
this approximation is acceptable. The results calculated 
by Eq. (fl"9]) are shown in Fig[TJ 

In comparison with the initial curve g^ c , in the su- 
perconducting curve g^ c there is a shift of maxima of 
the EPI towards lower energies by a value of the gap A. 
(Note that in the S — c — S contact for a transition into 
a superconducting state the position of maxima of the 
EPI does not change.) Besides, its amplitude is slightly 
smaller than that of the initial curve due to an additional 
broadening on a blurring factor S (Eq. (TTT10 . As it was 
already noted, in a superconducting state the spectrum 
of the EPI should show up in the first derivative of the ex- 
cess current. Indeed, if from the superconducting curve 
<7p C the initial normal curve g^ c is subtracted, then we 



4 




10 12 14 16 18 

eV, mV 



20 22 



FIG. 1: (a) gp C - point-contact function of EPI in Ta, re- 
constructed from the spectrum in Fig[2] <?p C - a theoretically 
predicted transformation of the point-contact function of EPI, 
upon a transition to a superconducting state (Eq. ([9]), see de- 
tails in the text); (b) is a difference curve; (c) is an integral 
from the difference curve, the PC function of EPI on the first 
derivative of an excess current. 



obtain the curve g^ c — g^ c (see FigQJb)). As it follows 
from Eq. (TT2"j) , the function of the EPI can be restored 
from the first derivative of the excess current, 
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This function is given by the curve shown in FigJIJc). 
Note that from Eq.© it follows that the point-contact 
function of EPI in the first derivative of the excess cur- 
rent is a superposition of g N (ui) and G s (u>). Thus, a 
position of phonon maxima in gp C (eV) is intermediate 
between g^ c and g!*, that is seen in Fig[TJ Besides the 
inelastic component of the superconducting addition in 
the spectrum, considered above, and taking only place in 
contacts with the noticeable excess current, for the point- 
contacts containing superconductors with a strong EPI 
an elastic addition to the spectrum is also possibl e) 10 ' 11 . 
This addition is due to an energy dependence of the su- 
perconducting order parameter A(e) and takes place for 
both point-contacts with direct conductivity (Eq.(lO) in 
Ref JlOT) and tunneling point-contacts (Eq.(ll) in RefilOl). 
In Refflll it is noted that for ballistic contacts nonlin- 



earity of the IVC of the elastic origin can be compara- 
ble to inelastic point-contact contributions. Obviously 
that for point-contacts with a large tunneling parameter 
Z it is exactly an elastic contribution will be prevalent. 
For both tunneling and point-contacts with direct con- 
ductivity the elastic addition to the spectrum does not 
contain the function of the EPI g(ui) in explicit form. 
However, it can be reconstructed by inverting the Eliash- 
berg equations (similar to the tunneling spectroscopy of 
Rowell-McMillan 12 ). For point-contacts with the direct 
conductivity or small value of Z it seems that the re- 
lation between elastic and inelastic contributions is de- 
termined by strength of the EPI. It is also necessary to 
take into account that for different superconductors as 
the EPI increases the energy mean-free path decreases 
leading to an increase in the concentration of nonequilib- 
rium phonons in the peripheral regions adjacent to the 
contact (see, for instance, RefQ). So for them one can ex- 
pect an increase in a relative part of nonlinearity, related 
to the inelastic superconducting scattering, in compar- 
ison with superconductors with a weak coupling. First 
experiments on measurements of point-contact spectra of 
Sn and In, satisfying the conditions ([1]) and @ a super- 
conducting state and their comparison with the spectra 
in a normal state have been performed in RefsQjHRj 
The authors have found excellent agreement between the 
Khlus theory^ and experimental results. At the same 
time experiments in Pb ( Ref 1171) have revealed a signif- 
icant contribution of elastic processes. By the way, in 
Fig. 2 in Ref j 1 71 one can see a form of the elastic contribu- 
tion to the spectrum, calculated by the abovementioned 
formula for the second derivative of the IVC of the Pb 
S — c — N point-contact, in comparison with experimental 
data. In the Sec. Ill we will restore a function of the EPI 
from the first derivative of the excess current assuming 
that the elastic addition can be neglected. The point- 
contact spectra of Ta, 2H-NbSe2 and MgB2 are used as 
objects of investigations. We estimate a contact diameter 
by the Veksler interpolation formula, 
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For the sake of convenience we list parameters of the 
objects under investigation: 

Ta (RefsQUU: - T C =4.47K; plra5.9-l{r 12 n-cm 2 ; 
£ =92nm; An=0.71meV. 

NbSe 2 (RefsiHg): - T C =7.2K; pl ab «2.2-lCT 11 r2-cm 2 ; 
£af>(0)=7.8nm; £ c (0)=2.6nm. Values of the supercon- 
ducting gaps within the two-band approximation (see, 
for instance, Eqs.(2) and (3) in Refl23l) estimated by 
us for the contacts given below with the resistance of 
1100 and 750 0, respectively: Ai=l.lmV, A 2 =2.46mV, 
k=0.62, A „ er =1.61mV; Ai=1.07mV, A 2 =2.48mV, 
k=0.7, A ot)er =1.49mV; Z=0, r 1>2 =0. 
MgB 2 (Ref El): T C =39K; p\ ab ^2A-W~ 12 n-cm 2 ; 

£ a6 (0)=12nm, £ c (0)=3.6nm, A i=1.8-=-2.4meV, 
A 02 =6.8-h7.4meV. 
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FIG. 2: (a)-Spectra of EPI of the Ta-Cu point contact in 
normal and superconducting states: T=1.6 K, Ro=70 fl, an 
initial dashed part of the superconducting curve containing 
the gap and nonequilibrium features is normalized by 100; 
(b)-a difference between superconducting and normal spec- 
tra, as well as an assumed form of a background curve; (c)-a 
difference curve (after background subtraction); (d)-a point- 
contact function of EPI reconstructed by integrating the pre- 
vious curve. 



III. EXPERIMENTAL DATA 
A. Point contacts Ta-Cu 

Fig|2] shows the PC spectra of the contacts Ta-Cu in 
normal and superconducting states. 

The resistivity of tantalum is p273A'=12.6-10~ 6 il-cm 
(Refill) the ratio of resistances pnzK I Pres ~18, which 
gives an elastic mean free path 84 nm in the electrode 
at the helium temperature. In this case the diameter of 
the point-contact with the resistance of ~70S1 is, accord- 
ing to Eq. ([21|). 3.8 nm. Note that in the estimation of 
the diameter the heterostructure was not taken into ac- 
count. However, considering that transport parameters 
of copper (pl«5.3-10 -12 f2-cm 2 (Refl26Tj) are close to those 
of tantalum, this correction can be neglected. A signif- 
icantly smaller background level in the superconducting 
spectrum is to be noted. An identical decrease of back- 
ground has been noted for the first time in RefQJ. In 
Refill it has been assumed that this effect can be related 
to decreasing the current of quasiparticles upon combin- 



ing some electrons responsible for the excess current into 
Cooper pairs during the Andreev reflection, since the EPI 
in the contact is connected with a quasiparticle compo- 
nent of the total current. It is necessary to note here 
that the shown spectrum satisfies the conditions (fTJ) and 
([2]) and is similar to most intensive spectra described in 
Nevertheless, upon a transition to a superconduct- 
ing state noticeable distinctions from predictions of the 
theory are observed. Here we will not discuss an appear- 
ance of nonequilibrium characteristics in spectra, which 
have been considered in detail in RefsHil and restrict 
our attention to the EPI. First of all we note that the 
greatest distinctions are observed in the range of small 
energies. A soft mode in the vicinity of 7-^8 mV starts to 
show up as a peak. Instead of broadening, the peak close 
to 11-^12 mV becomes sharper abruptly, but the distinc- 
tions in the vicinity of 18 mV are already not so pro- 
nounced. The sharpening of the peaks can be caused by a 
high concentration of nonequilibrium phonons with small 
group velocities which leave the contact region slowly. 
Their velocities correspond to the frequencies at which 
the density of states is maximal. The probability of reab- 
sorption of phonons by electrons exhibiting the Andreev 
reflection increases with the phonons concentration. It 
is just a selection of phonons according to duj/dq = 
determines the effect of stabilization of peaks positions 
and their sharpening upon a transition into a supercon- 
ducting stated. In an N state phonons participate in 
inelastic backscattering processes in a broader area of 
the phase space, and the position of maxima in the spec- 
trum can change depending on orientation of the contact 
axis with respect to the crystallographic axes. In addi- 
tion, the greatest distinction of the curves in Figs. FigQ] 
and Fig [2] exactly in the region of small displacements 
can be due to the fact that an effect of the superconduct- 
ing form- factor entering to G (to) was not taken into ac- 
count in calculating the convolution in Fig|T] (curve <?p C ), 
see Eq. I|7|). As it was mentioned in Introduction it is 
just in the region of small displacements on the contact 
an effect of the superconducting form-factor manifests 
itself most pronouncedly, and the main contribution to 
the current is made by processes of scattering without 
changing the z-component of the velocity, of the type of 
the Andreev reflection. The absence of the sharpening 
of phonon peaks in the PC spectra of In and Sn upon a 
transition to a superconducting stated— is an indirect 
indication of a small contribution of the region close to 
the contact in the inelastic superconducting spectrum. 
For these superconductors the energy relaxation length 
at energies of phonon peaks is higher than that in Ta, 
and in spite of a small group velocity of such phonons 
their concentration in edges in the vicinity of contrac- 
tions is insignificant to affect noticeably a form of the 
spectrum. At the same time in the region of intensive 
generation of phonons the selection on the basis of the 
phonon group velocity is not so efficient due to a small 
size of the current concentration region. The presence 
in the S spectrum of the gap and nonequilibrium fea- 



6 



tures as well as changing the form of the spectrum and 
the reduced background makes it not possible to restore 



the function of the EPI from the curve V? 



JY 



simple integration. So, to compensate the influence of 
all neglected factors before integrating it is necessary to 
subtract the background as it is done for restoring the 
function of the EPI from a spectrum of a normal state. 
Since at energies above the Debye ones the function of 
the EPI turns to zero, the curve after subtracting the 
background should become to be zero at these en- 
ergies. The second necessary condition is an absence of 
a background on the curve obtained as a result of inte- 
gration. To fulfill this condition the areas under the 
curve after background subtraction should be the 
same above and below an abscissa axis. Note that 
above mentioned requirements are satisfied for different 
forms of a background curve. However, the variations do 
not lead to significant changes of the form and positions 
of phonon characteristics of a restored function of the 
EPI gp C pc. So for a rough estimate of the EPI func- 
tion, it is enough to plot a background curve freehand. 
Because the gap and nonequilibrium features do not al- 
low to reconstruct correctly a trend of the curve in the 
range from to ~6 mV, this part was approximated by 
a parabola segment. In order to prevent from additional 
misrepresentations the background curve was smoothed 
out by segments of polynomial curves. We note in pass- 
ing that the restored function of the EPI gp C is close 
enough in form to the EPI function of a superconducting 
state, which could be obtained from the ^-spectrum in 
FigJ2] (a) by subtracting a traditional background. 



B. Point-contacts 2H-NbSe 2 -Cu 

The investigation of the EPI spectra of 2H-NbSe2 in a 
superconducting state is presented in Refl3ll. We empha- 
size here that this compound is very difficult for exper- 
iments due to a strong dependence of the point-contact 
spectra on the quality of initial samples of single crys- 
tals. For example, one of the first single crystals investi- 
gated by us had a resistance ratio P3oo/p r es ^30. Point- 
contacts based on these single crystals were character- 
ized by a very low noise level up to the resistance of 1 
kil. At the same time the point-contacts based on single 
crystals with P3oo/Pres ~100 and having the resistance 
higher than 10-^20f2, in most cases had a very high level 
of noise that did not allow to use them for studying the 
EPI. A possible reason for this is incommensurate charge 
density waves (CDWs) that exist in this compound. In 
the former case the single crystals are not perfect enough, 
and there are effective pinning centers, preventing their 
movement. In the later case there are not many of such 
centers, and the motion of the CDWs in the vicinity of the 
contact in the current state generates such noise. Note a 
principal difference between the point-contact EPI spec- 
tra in 2H-NbSe 2 FigEJand Ta FigH 

In a normal state there is no structure in the region 




FIG. 3: (a) and (b) - Point-contact spectra of 2H-NbSe 2 -Cu 
and assumed background curves, T=4.2 K, H=0. The contact 
with the resistance of 750S1 is obtained from the previous 
one by shortening, (c) and (d) - Spectra after background 
subtraction 



of phonon energies in the spectra of 2H-NbSe2 (see Fig.l 
in Ref.3I). At the same time in a superconducting state 
the spectra have many features. It is noteworthy that by 
the standards of the traditional PC spectroscopy contacts 
are very high-ohmic. The reasons for selecting such con- 
tacts are a high resistivity of the compound, greater by 
about 2 orders of magnitude than in conventional metals 
Oil ~2-10- 4 fi-cm; p x ~l(r 3 fi-cm (ReiH)), and, ap- 
parently, a small energy mean free path of electrons at 
energies close to the Debye due to the strong EPI. The 
reason for the phonon features to appear only in a super- 
conducting state is apparently associated with a tech- 
nology of manufacturing the PCs. The sample of 2H- 
NbSe2 was a single-crystalline plate with the dimensions 
of ^15x20x0.1 mm. The electrodes with the dimen- 
sions ~3x5 mm were cut out by a blade and glued with 
silver paste to a wire holder. Thus, the contact pads 
of the electrodes were machined, and a crystal struc- 
ture of 2H-NbSe2 on a site of the cut was largely bro- 
ken. Consequently, the point-contact is very inhomo- 
geneous; there are clean edges and a dirty core with a 
strongly distorted lattice with a short mean free path of 
electrons. In a normal state, an intensity of the phonon 
spectrum is proportional to the ratio k/d and very small. 
The spectrum itself is formed by a material in a volume 
of the order of a diffusion mean free path A 6 and, be- 
cause of the crystal structure distortion, contains only 
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a background component. An estimate of diameters of 
the contacts by Eq. (j2"Tj) gives 2-^3.4 nm for the first and 
2.3-^4.3 nm for the second contacts respectively. An elas- 
tic mean free path for these estimates was about 33 and 
1.1 nm, corresponding in the former case to the resistiv- 
ity of unperturbed material in the vicinity of Tc and to 
p in 2H-NbSe2 at room temperature in the latter case. 
In the other words, it was assumed that in the contact 
region p^ooPres Taking into account the lack of spec- 
tra in a normal state, we can assume that the mean free 
path near the narrowing is even smaller, and diameters 
of the contacts are larger. Upon a transition to a super- 
conducting state at the IVC of the point-contact there 
appears the excess current which is formed by processes 
of the Andreev reflection at distances of the order of the 
coherence length. This length is larger than the diame- 
ter of the contact, and the greatest contribution to the 
excess current is given by those of the edges where the 
crystal lattice is the most perfect, because the defects 
make worse superconducting parameters. Under this as- 
sumption there is only dirty area of the superconducting 
electrode, and the second normal electrode remains clean 
in this case. This fact provides a unique opportunity to 
obtain spectra that reflect a fine structure of the EPI of 
an undistorted crystal lattice. As noted above, the su- 
perconducting PC function of the EPI is proportional to 
a function of the first derivative of the excess current. 
Restoring the function of the EPI from a superconduct- 
ing spectrum in the absence of any features on the second 
derivative in a normal state assumes only subtraction of 
a monotonic background curve. In this case, in the vicin- 
ity of the gap region (~3-^4 mV) the background curves 
were plotted so that initial parts of the difference curves 
(FigOl curves c, d) were parabolic. The point-contact 
functions of the EPI restored from the spectra c and d 
(Figj3|) via integration are shown in Fig|4] 

Note that the available literature data on the Debye 
energy of a phonon spectrum vary greatly among them- 
selves. The only known calculated function of the phonon 
density of stateS ) 33 i 34 according to the authors, is based 
on elastic moduli and the maximum frequency of the 
quasiparticle mode and is defined up to the energy of 
^32 meV. In this interval, it agrees well with our data. 
One of the co-authorS) 33 i 34 S.B. Feodos'ev, suggested to 
associate the features observed at higher energies with 
two-phonon processes, the calculated function of which 
is shown in the same figure. The convolution presented 
coincides perfectly in a form with the data obtained by 
us, however so high intensity of two-phonon processes is 
unclear to us. On the other hand, for much higher fre- 
quencies the peaks observed on our curves are in good 
agreement with flat parts of the calculated dispersion 
curves, where <9w/<9q=0 in the Feldman theoretical model 
for a dynamic matrix 35 (in the figure the first model is 
presented), as well as with the calculations in Refl36l. 
The exception is the most high-energy peak at ~60 mV. 
Possible causes of this peak are discussed in Refl3ll A 
position of the phonon features in our spectra correlates 




FIG. 4: A comparison of point-contact functions of EPI for 
2H-NbSe2, reconstructed from the spectra (c) and (d) in Fig[3] 
3, with a function of the phonon density of states from Refsl3a 
and Hi] (dash) as well as with a twophonon convolution (dot) 
(the lower part of the figure) and the calculated dispersion 
curves (model 1) from Refl35l (for details see the text). 



well with the data of experimental works on the inelas- 
tic neutron scattering and the Raman spectra. Unfortu- 
nately, the experimental neutron data for this compound 
is generally limited to the study of low-energy phonon 
branches. The small variations observed in the position 
of the phonon features at ~l^-2 mV are apparently re- 
lated both to an accuracy of the theoretical model and 
an influence of the superconducting energy gap in our 
data. Note that recently the point-contact functions of 
the EPI have been obtained for the 2H-TaSe2-compound 
belonging to the same crystallographic group 3 -^. It is in- 
teresting that a boundary of the phonon spectrum in this 
compound is greater than 70 mV. 



C. Point-contacts MgB2-Ag 

For the analysis, the superconducting PC spectra of 
MgB 2 (Fig©, reported in RefslMffl were used. 

The spectra were obtained by touching a massive sil- 
ver electrode to the film of MgB2 and represent an inter- 
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FIG. 5: Point-contact spectra of MgB 2 (T=4.2 K, H=0), 
taken from Refsl38M4"ol. as well as assumed background curves. 
For clarity, the curves are shifted vertically. 
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FIG. 6: Point-contact spectra of MgB2 (Fig[5]) after subtract- 
ing background. Parts of the curve up to ~18 mV are ap- 
proximated by a parabola. 





mediate case between Ta and 2H-NbSe2 for the appear- 
ance of features in a normal state. For the point-contacts 
similar to those shown in FigJSJ in the normal state the 
weakly pronounced highly blurred features are observed 
in the region of the phonon energies with a high-level 
background^. The residual resistivity of the film was 
6T0 _6 f2-cm, that corresponds to the elastic mean free 
path of 4 nm. The diameters of the given contacts are 
estimated to be 4.4, 3.9, 3.8, 2.2, and 1.9 nm (36, 43, 45, 
110, 150 f2) respectively. Thus, the estimates are very 
close to the estimated diameter of 2H-NbSe2, and the 
values of the elastic mean free path in the contact region 
in the MgB2 is found to be several times larger. Such es- 
timates are quite reasonable, considering that we almost 
do not make additional lattice distortions. When a back- 
ground is subtracted form the curves in a superconduct- 
ing state, we were guided by the conditions specified in 
SubSec. A and highlighted by bold. FigJS] shows the sec- 
ond derivatives of the IVC and the background curves, 
and FigJB] demonstrates the difference curves. 

As can be seen, the background curves are not mono- 
tonic as in the case of 2H-NbSe2, and contain one or 
two blurred maxima. Note that presented in Reffill 
point-contact spectra of the EPI in a normal state as 
a whole are similar to the curves represented by the 



background curves and contain the blurred maximum at 
around ~60 mV. The point-contact spectra of M2B2 in 
a superconducting state were interpreted in Refj4fj from 
the point of view of an energy dependence of the super- 
conducting order parameter A(e). The EPI function re- 
constructed from the difference curves are shown in FigJT] 

Fig[5] demonstrates the PC function of the EPI av- 
eraged over five contact s, g p™ r , related to it the calcu- 
lated g ca ; c (eV) (Refsl42|-|44l) and obtained from the tun- 
neling data g ea; p(eV) (RefsH and l45h functions of the 
EPI as well as the functions of the phonon density of 
states F exp (eV) and F co { c (eV)i§ 

Unfortunately, it is impossible to compare in detail 
this function with the functions reconstructed from the 
PC spectra in a normal stated due to a rather diffuse 
structure of the latter. So we compare with the point- 
contact function gpc(eV) for ZrP>2 (Refl47l) belonging to 
the same crystallographic group. As can be seen, the 
functions are very similar in shape and position of the 
phonon features which for the averaged spectrum are at 
the energies by 4-^5 mV smaller than for ZrB-2, with the 
exception of the first peak shifted to the higher ener- 
gies by ~16 mV. This is understandable, since this mode 
is associated with vibrations of the Zr ions, the atomic 
mass of which is ~3.75 times greater than that of the 
Mg ions. The calculated functions of the EPI g ca ic(eV) 
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FIG. 7: Point-contact functions of EPI in MgE>2, recon- 
structed from the difference curves in Fig(6] 

(Refs fi^ - li^ are notably different in form. However, their 
comparison with the average PC function gp™ r demon- 
strates a close value in the energy of individual branches 
of the spectrum. As for the tunneling function of the 
EPI geap(eV) ) 44 ' 45 it can be assumed that the maximum 
at 58 mV corresponds to the peaks at 46 and 66 mV of 
our averaged functions. These peaks for the tunneling 
contact could not be resolved due to insufficient quality 
of the superconducting film. This assumption can also 
be supported by softening the high-energy peaks from 90 
to 85 mV. A comparison of the function of the phonon 
density of states, both calculated F CQ / c (eV) and obtained 
from the neutron data ¥ exp (eY)^ with our demonstrates 
a good correlation in a position of phonon peaks. In 
general, it can be concluded that the method used can 
accurately restore a form of the EPI function from su- 
perconducting point-contact spectra. 




10 20 30 40 50 60 70 80 90 100110 
eV, meV 

FIG. 8: A comparison of the point-contact function of EPI in 
MgB2, gpc r ( a ), averaged over five spectra (FigO: with cal- 
culated functions of EPI, g ca !c(eV) fRefs l43 - l44[ ) (b); functions 
of the phonon density of states f Ref l46h . calculated F ca ; c (eV) 
and experimental F exp (eV), obtained from neutron experi- 
ments (c); the PC function of EPIgpc in ZrB 2 (Refill) (d); 
the EPI function g exp (eV) (RefsO and l45h . reconstructed 
from tunneling experiments (e). 



an interpretation of superconducting spectra is most 
straightforward in two limiting cases: 

1. for ballistic point-contacts in which the relative 
contribution to the spectrum of a superconducting 
component, as follows from the Khlus theory— is 
small and retains the overall spectrum shape, suf- 
ficiently close to a normal state; 

2. for point-contacts with a very dirty core and clean 
edges. 



IV. DISCUSSION OF THE RESULTS 

As experiments demonstrate the method is effective 
for metals both with a large coherence length and with 
rather short one. For clarity, the characteristic param- 
eters of the point-contacts considered in this paper are 
given in the summary table. 

Based on the above material, we can conclude that 



Theoretically, a decrease of the elastic mean free path 
in the contact region should lead only to a decrease in 
the intensity of the spectrum and should not affect its 
shape. But in reality, such a decrease of the mean free 
path is caused by either impurities or defects and distor- 
tions of a crystal lattice. In any case, this is a violation of 
short-range order, accompanied by a blurring the spec- 
trum. In the limiting case for amorphization of a contact 
region the spectrum may be absent completely, but the 
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TABLE I: Parameters of point-contacts Ta, NbSe 2 , MgB 2 . 





R, Q 


pi, fi • cm 2 


P30a/pres 


pres, f2 • CHI 


d, nm 


k, nm 




A, mV 


Ta 


70 


5.9-10" 12 


18 


0.7-10" 6 


3.8 


84 


92 


0.71 


NbSe 2 


750-M000 


2.2-10" 11 


30 


6.7-10" 6 


2H-2.3 


33 


7.7 


1.1^-2.5 








1 


2-10" 4 


3.4^-4.3 


1.1 






MgB 2 


36^150 


2.4-10" 12 


2.2 


6-10~ 6 


2.2^4.4 


4 


12 


2^7 



contact will not go into a heating regime. If edges at the 
same time remain clean, after subtracting a background 
the spectrum is easily obtained from the first derivative 
of the excess current. Note that to realize the situation 
of dirty narrowing and clean edges is quite difficult. Most 
often, edges also appear to be with the short-range order 
breaking. A theory for the diffusion limit is presented in 
Refl48l. The solution is obtained only for a long channel 
model (bridge), provided that the bridge and the adja- 
cent edge areas contain large amounts of impurities. In 
the diffusion limit, as in pure point-contacts, there are 
processes at which electrons injected into a superconduc- 
tor in the process of the Andreev reflection interact with 
nonequilibrium phonons, giving rise to a "back current." 
In addition, the contribution to the conductivity is made 
by processes such as relaxation of charge imbalance and 
the virtual processes that are sensitive to the frequency 
dependence of the EPI function. In Refl48l it is empha- 
sized that these additional channels are observed only 
for the geometry of a dirty bridge with a massive super- 
conductor. Obviously, in most cases, this geometry has 
little to do with a real experimental situation. Depart- 
ing from the bridge geometry and considering the actual 
spatial dependence of all quantities in the problem for 
a transition from a narrowing to massive edges compli- 
cates significantly the problem and still has not found 
a solution. Nevertheless, it is obvious that a decrease 
of the elastic mean free path in the vicinity of a nar- 
rowing reduces much faster a normal component of the 
spectrum than superconducting. Therefore, upon a tran- 
sition from normal to superconducting state an intensity 
of the spectra of the EPI increases in such point-contacts. 
An example of such a transformation of spectra can be 
observed in Refs l49ll50l . An intermediate regime, close to 
diffusion one, is very difficult to interpret. In this regime, 
both contributions are similar in intensity, and the nor- 
mal spectrum has a small blur. A superposition of contri- 
butions of similar intensity greatly distorts a shape of the 
resulting spectrum, and separation of its superconduct- 
ing part without carrying the contact to a normal state 
is a highly non-trivial task^i. A complicating factor is a 
possibility of appearing on the derivatives of the CVC the 
nonspectral features of the thermal and nonequilibrium 
origin, related to the destruction of superconductivity in 
the contact region due to the spatial inhomogeneity of a 
superconductor. To exclude from the consideration such 
features there should be some statistics of the spectra, 
since these features are not sufficiently reproducible. In 



addition, in contrast to the spectral, their position on 
an energy axis, are generally more dependent on tem- 
perature. A position of the spectral features of the ex- 
cess current on an energy axis depends on the gap value 
and, hence, the temperature. However, taking into ac- 
count that the gap starts to decrease significantly only 
at T>2/3Tc, in some temperature range this position 
can be considered practically constant. In addition to 
the above-noted difficulty of separating the normal and 
superconducting contributions to the spectrum and ex- 
cluding the nonspectral features, the empirical method 
of subtracting the superconducting background is rather 
arbitrary yet. The problem of obtaining quantitative in- 
formation from the spectrum, especially in the case of 
irregular distribution of impurities, is not solved. While 
in the simplest cases for materials with a large coher- 
ence length, the purity of the contact can be judged by a 
value of the excess current, in the compounds with small 
£ the excess current is influenced by other mechanisms 
in addition to the Andreev reflection. Therefore, in com- 
plicated cases of the greatest practical interest, such as 
high-temperature superconductors, this information will 
only be of a provisional character. Note that in the case 
of a short coherence length and a relatively low critical 
current density in the current state, in the vicinity of the 
point-contact there can be phase-slip surfaces preventing 
the observation of phonon features^. Thus, though the 
PC superconducting spectroscopy is not as universal as 
a traditional, in some cases it is more efficient. First of 
all, these are compounds which are difficult to take to a 
normal state, as well as compounds in which in a normal 
state the spectra of the EPI are generally of poor quality 
or absent at all, as, for example, in NbSe2- 



V. BRIEF SUMMARY 

1. Using a Ta point-contact as an example it was 
shown that it is possible in principle to reconstruct 
the point-conact function of the EPI from the first 
derivative of the excess current. It was found that 
the restored PC function of the EPI differs in form 
from the EPI function of a normal state, and the 
largest differences are observed at low energies. An 
explanation of the observed differences was pro- 
posed. 

2. From the PC spectra in a superconducting state, 
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for the first time the function of the EPI was re- 
constructed in 2H-NbSe 2 . A position of phonon 
peaks in the EPI function agrees well with the flat 
parts of calculated dispersion curves. 

3. The PC functions of the EPI in MgB 2 , restored 
from superconducting spectra for the first time, 
show good reproducibility. A comparison with re- 
lated functions of the density of phonon states, 
F(x), tunneling and calculated functions of the EPI 



shows a small shift of main phonon peaks to lower 
energies in full agreement with theoretical predic- 
tions. 

In conclusion, we wish to thank Yu. G. Naidyuk for 
numerous tips, comments and additions made during dis- 
cussion of the paper. Our thanks are also to S. B. Feo- 
dos'ev for a discussion of the phonon spectrum of 2H- 
NbSe 2 . 
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